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Second order perturbations in the radius stabilized 
Randall-Sundrum two branes model 
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The nonlinear gravitational interaction is investigated in the Randall-Sundrum two branes model 
with the radius stabilization mechanism. As the stabilization model, we assume a single scalar field 
that has a potential in the bulk and a potential on each brane. We develop a formulation of the 
second order gravitational perturbations under the assumption of a static and axial-symmetric five¬ 
dimensional metric that is spherically symmetric in the four-dimensional sense. After deriving the 
formal solutions for the perturbations, we discuss the gravity on each brane induced by the matter 
on its own side, taking the limit of large coupling of the scalar field interaction term on the branes. 
We show using the Goldberger-Wise stabilization model that four-dimensional Einstein gravity is 
approximately recovered in the second order perturbations. 
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I. INTRODUCTION 

Many unified theories require spacetime dimensions higher than that observed in the Universe, and thus the extra 
dimensions must be invisible by some mechanism. One of the possible schemes is known as Kaluza-Klein compacti- 
fication. Recently, theories with extra dimensions have attracted considerable attention from the other viewpoint of 
providing a solution to the hierarchy problem [1-3] [4,5]. The main idea to resolve the large hierarchy is that the small 
coupling of four-dimensional gravity is generated by the large physical volume of extra dimensions. These theories 
provide a novel setting for discussing phenomenological, cosmological, and conceptual issues that are related to extra 
dimensions. 

The model that was introduced by Randall and Sundrum (RS) is particularly attractive. The RS two branes 
model is constructed in a five-dimensional anti-de Sitter (AdS) spacetime [2]. The fifth coordinate is compactified on 

jZi, and the positive and negative tension branes are on the two fixed points. It is assumed that all matter fields 
are confined on each brane and only the gravity propagates freely in the five-dimensional bulk. In this model, the 
hierarchy problem is resolved on the brane with negative tension if the separation of the branes is about 37 times the 
AdS radius. 

Apart from the fine tuning of the brane tension that is necessary to solve the cosmological constant problem, one 
of the significant points in discussing the consistency of this model is whether four-dimensional Einstein gravity is 
recovered on the brane [6-12]. Another point is to give a so-called radius stabilization mechanism that works to 
select the required separation distance between the two branes to resolve the hierarchy without fine-tuning [13-23]. 
The stabilization mechanism not only is important to guarantee a stable hierarchy, but also plays an important role 
in the recovery of four-dimensional Einstein gravity in linear order [7,12] and of the correct cosmological expansion 
[24-28]. The discussion of the recovery is almost independent of the details of the stabilization model. The essence in 
recovering the four-dimensional Einstein gravity in linear perturbations is that the massless mode of the scalar-type 
gravitational perturbation disappears due to the bulk scalar field, and only the tensor-type perturbation continues to 
have a massless mode. When the stabilization mechanism is turned off, the induced gravity on the brane becomes of 
the Brans-Dicke type with an unacceptable Brans-Dicke parameter [6]. 

A large number of studies have been made of the gravity in the brane world model [29-37]. Although the model 
does not have a drawback in linear perturbation, it is not a trivial question whether the second order gravitational 
perturbation works as well. For the second order perturbation in the RS single brane model without bulk scalar 
field, where the tension of the brane is positive, it has been confirmed that there is no observable disagreement 
with four-dimensional Einstein gravity [38,39]. However the setting of the RS two brane model with stabilization 
mechanism is quite different from that of the single brane model, and furthermore we are mainly concerned with the 
gravity on the negative tension brane. In this paper we study the second order gravitational perturbation of the RS 
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two branes model with a stabilization mechanism due to a bulk scalar field. To simplify the analysis, we consider 
static and axisymmetric configurations, which means that the metric on the branes is spherically symmetric. After 
developing a formulation to calculate the second order perturbation, we take the limit that the coupling of the scalar 
field interaction term on each brane is very large. In this limit, we find that four-dimensional Einstein gravity is 
approximately recovered. 

The paper is organized as follows. In Sec. II we describe the model that we will study, and derive the second order 
perturbation equations in the five-dimensional bulk. We also discuss the gauge transformations and the boundary 
conditions. In Sec. Ill we explain our approximation scheme, and give the formal solutions. In Sec. IV we review the 
results for linear perturbations, giving their explicit expressions following the notation of this paper, and explain the 
setup of the problem that we study in the present paper. In Sec. V we analyze the second order metric perturbations 
induced on each brane. We show that the four-dimensional Einstein gravity is recovered with some small corrections. 
These results are summarized in Sec. VI. 


II. PERTURBATION EQUATIONS IN THE RS MODEL 

We consider the second order perturbations in the RS two branes model with a five-dimensional scalar field intro¬ 
duced to stabilize the distance between the two branes. According to the warped compactification of the RS model, 
the unperturbed metric is supposed to be 

ds^ = gabdx°'dx^ = dy"^ + a^{y)r]f_i,ydx^dx ‘', ( 2 . 1 ) 

where rj^i, is the four-dimensional Minkowski metric with (— -I- -I- -f) signature. The y direction is bounded by two 
branes located at y = y+ and y = y-, whose tensions are assumed to be positive (A(_|_) > 0) and negative (A(_) < 0), 
respectively. On these two branes, Z 2 symmetry is imposed, and we adopt the convention y+ < y-. To generate the 
hierarchy between Planck and electroweak scales, we need 

( 2 . 2 ) 

a_ 


where a± = a{y±). 

In this paper we investigate the gravity induced by non-relativistic matter fields confined on each brane whose 
energy-momentum tensor is given in the perfect fluid form: 

= ai4diag{-p±,P±,P±,P±}. (2.3) 

The warp factor in the definition of the energy-momentum tensor (2.3) is incorporated for the following reason. In 
the present analysis, we adopt the normalization that any physical quantities are always mapped onto and measured 
by the length scale at y = y+. Since the length scale is warped by a warp factor a{y), physical quantities such as p_ 
and P- are scaled by a factor al^. 

To simplify the analysis, we restrict our consideration to a static and axisymmetric spacetime whose axis of symmetry 
lies along y direction. We denote the perturbed metric by gab = gab + hab- The four-dimensional perturbation hfj. 1 , is 
divided into a trace part and a transverse-traceless (TT) part. According to this decomposition, we assume that the 
perturbed metric has the diagonal form 


ds^ = -I- a'^{y) 


_^Mi~yV)-‘>P(.r,y) ^^2 _|_ ^B{r,y )-tp{r,y) ^^2 _|_ ^C{r,y)-'ip(r,y) ^2 ^^2 


(2.4) 


Here A, B, and C correspond to the TT part, and ^ to the trace part. The TT condition to the linear order is given 
in terms of A, B, and C as 


A^'^\r,y) = -^dr{r^B'^^\r,y)) , 

C'«(r,2/) = l5.(r2pW(r,y)), (2.5) 

2r 

where we have expanded the metric functions up to second order as 

A= (2.6) 

J^l,2 
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The other metric functions V and ip are expanded in the same way. Henceforth we neglect higher order terms without 
mentioning it, and we omit the superscript indicating the order when it is obvious. We impose the same conditions 
as Eq. (2.5) on and so that and are derived from once it is solved. Since the trace of 

the metric is given by 

r"Ku = - 2V>" - + 52 + 2C^), (2.7) 

the trace part at linear order is while does not correspond to the trace part at second order. Hence the 
second order counterpart of the condition (2.5) does not mean that is the transverse-traceless perturbation, but 
we extensively refer to these metric functions as the TT part and ip as the trace part. Later we will show that Y 
coincides with p) 'm & linear perturbation, and thus this metric assumption is the same at least in linear perturbation 
as the “Newton gauge” condition. 

The Lagrangian for the bulk scalar field is 

- Vb((^) - ^ V(„){0)5{y - Va) , (2.8) 

a=± 

where Vb and V(±) are the potential in the bulk and that on the corresponding brane. For most of the present analysis, 
we do not need to specify the explicit form of the potentials Vb and V(±). The scalar field is expanded up to second 
order as 

(p{r, y) = (po {y) + (r, y) + (r, y ), (2.9) 

where (pQ is the background scalar field configuration, which depends only on y. 

From the five-dimensional Einstein equations with the cosmological term A and the equation of motion for the 
scalar field, we obtain the background equations as 

Hiy) = -| 0 g(j/), 

f (^<^ 0 ( 2 /) - VBiMy)) - 

^o{y) + 4.H{y)(poiy) - V^(0o(2/)) = 0, (2.10) 

where 7J(t/) := a(j/)/a(j/)(« —a/—A/ 6 ) and the five-dimensional Newton’s constant is G 5 = k/8tt. An overdot denotes 
differentiation with respect to y. 

A. 5D Einstein equations in the bulk 

In this subsection we derive the master equations for and from the five-dimensional Einstein equations. 


From the (r, r), {9,9), and {r,y) components of the Einstein equations, we obtain two independent equations: 

= -^{dyi;^-^^ + 2HY^^'> + e^^'>Sy,), ( 2 . 11 ) 

2k(Po 

where A = Here and Sy, are the second order source terms that are constructed from the first order 

quantities, and the explicit forms are given below. We have introduced the symbol that is defined by = 0 and 
e(2) = X to represent the first and the second order equations in a single expression. These equations are reduced to 

^(■^)(r,2/) (2.12) 

(p('^)(r, y) = ^ dy{a^Y^-’'’) + —-f dyA-^sJ . (2.13) 

2/s:0o<^ 2/’C0o ^ ^ 

The equations for A^'^^ and yGl are obtained from the (r, r) and {y,y) components of Einstein equations as 

r^CTT) 1 A](a2A(‘^)) = , (2.14) 

[L^^'> + , (2.15) 
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where 


:= ^dya^dy^ , 

L(’^) := a^^ldy^^dyO^ - . 

a^n 


The second order source terms Sa and Sy are given later. After the simplification using the linear order equations 
(2.12), (2.13), (2.14), and (2.15), 5^ and are written down as 


drS^{r, y) = ^dr(r^(ABf) + + ^rdrAB - 


H- 2^y * ^.r^.v H“ 


r-^ ,y 


2^8/3 




and 


where 


= — 




a‘^4'0 ^ 


■ (ay- 


Ka'^(j)Q 


-ip 


(2.16) 


(2.17) 


S^ = 


^ J (^B^yA^r + ^B^yr{3A^r + B^ri^dr — J B^yY^rdr -\—J p^rAYdr + HY'^—^dy[p^). (2-18) 


In the derivation of Eq. (2.16), we used Eq. (B14). The underline in Eq. (2.16) is attached for convenience of our 
explanation. So is the double underline below. Using Eq. (2.16), Sa is given by 


5a = --BB,r - ^{B^rf - ^dr{r’^BB,rr) " /dr(^(B„)2 + -(B,)^) 
r o r’^ J \2 r J 

driYAAYlA , ' j driaXrBA-2 j 


-3 


r V Ka'^(j)Q 


-,2'^V 


4>o 


(2.19) 


where we again used Eq. (B14). The complete expression for Sy is slightly complicated. For later convenience we 
divide Sy as 


Sy = Sy — (j)Qa^ dy 


1 r 


+ dyA ^ S'tp 


where §y is given by 

,2 


§y = -^((A.,)2 + (B,)2 + ^A,B,,) - y'|AB(yAi? + ^(AB),, + 4i?,,) +^i?,,(Ai?),,+ 


iB,r 


3r 


a^H 


r / X \ 

A^r-{3A + B)^y-rB^rB^ry + Yr(AB+ — d:\-B{AY)^. 

+ (y + {AY)^r + - p{Ap)^r) “ Idr 

V ())g / 3 \ Av) 


-2BrAY 


9(Ay)- 


+ 3Y^rAY 
- a^HY^ . 


( 2 . 20 ) 


( 2 . 21 ) 


B. Boundary condition 

In the previous subsection, we derived master equations for the metric functions in the bulk up to second order. To 
solve these equations we must determine the boundary conditions on the branes. It is well known that the boundary 
condition is given by Israel’s junction condition [40], which is easily obtained in Gaussian normal coordinates. On the 
other hand, the Newton gauge simplifies the master equations for the perturbations. Therefore we consider the gauge 
transformation between them. 

In Gaussian normal coordinates the metric becomes 
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ds"^ = dy^ + a^{y) -e^dP + e^dP + e^PdVP 


(2.22) 


with y = const on either the positive or negative tension brane. Note that we introduced two sets of Gaussian 
normal coordinates; one is the coordinate set in which the positive tension brane is located at y = y+, and the other 
is that in which the negative tension brane is located a,t y = y-. Corresponding to these two different Gaussian 
normal coordinates, there are two infinitesimal gauge transformations x°' = x°' + between the Newton gauge 

and the Gaussian normal gauge, respectively. To satisfy the restriction on the metric form in both gauges, the gauge 
( 1 ) ( 2 ) 

parameters associated with each brane must take the form of 


(J) py / 1 r 1 ( 1 ) 1 \ 

(r, y)=J dy' j + (r), 

(J) fy f 1 c-') PJ) ( 1 ) r- (1) ‘'P 

(r, y) = dy'i -- e±,r +— e±,r B-Y + 2He± + C±,r + C± {r), 
d y± \ ^ ® 


(2.23) 


(2.24) 


where we simplified the integrand of the equation for the second order perturbation by using the result for a linear 
perturbation. The functions of r, ^5, arise as integration constants. The arbitrariness of is due to a residual 

gauge degree of freedom of the coordinate transformation in the radial direction, while is determined with the aid 
of the junction conditions as we will see below. 

The gauge transformations for each metric component are given by 

(J) |-(1) _ (1) _ (1) n 

A^^\r,y) = - 2H , 

(J) r(i) _ (1) _ (1) 1 1 

B^^Hr,y) = - 2H e± -2 5^ + Ci ^ , 


(J) 9 (J) r^ij (i; I (1) (,i; 1 

C^^\r,y) = C'(^) - - 2H e± --C± G + Hie±? 


(2.25) 


As for the scalar field, its gauge transformation is given by 


My) = My), 

(r, y) = M'’ {r, y) - SM'’ (a y), 


where 


/ f-i-J (i) . 1 r . (1) 

SMHr,y) = <i>o K± - <^o ^±).r - -[0o(^±)^ 


(2.26) 


(2.27) 


As mentioned earlier, we assume the energy-momentum tensor to be of the perfect fluid form (2.3). The four¬ 
dimensional energy-momentum conservation = 0 becomes 


{p± + P±)drA^'^\r,y±)+ 2drP± =0, 


(2.28) 


and hence we find that P± is a second order quantity. This equation represents the force balance between pressure 
and gravity acting on the matter field. 

Now we consider the boundary conditions. Israel’s junction conditions on the 3-branes are given by 


= -^[M - ^ (at y = y±).. 


(2.29) 


where A(±) is the tension on each brane, and is the energy-momentum tensor for the scalar field. Here and 
hereafter, when we evaluate the value a,t y = y±, we take the value at j/ = j/± ± e in the e ^ 0 limit. Since, by 
assumption, the scalar field dose not have a kinetic term on the brane, its energy-momentum tensor on the brane is 
given by = —V(±)5)). Then, Eq. (2.29) gives at the lowest order 


TH = |y(±) + |A(±) {y = y±). 


(2.30) 
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The potential V(±) must be chosen to satisfy this condition. The (t,t) component of the junction condition (2.29) 
gives 

±Q^A^J) = = y^^ . ( 2 . 31 ) 

The junction condition in the Newton gauge is obtained by applying the gauge transformation to this equation. We 
define as the jump of the derivative of a metric function in the Newton gauge: 




y ly=v±’ 

l^\A) 


= -«(T°--r)^ iy = y±).. 


(2.32) 


where is given by the substitution of Eqs. (2.10), (2.15), (2.17), and (2.26) as 

2/^ / -7 1 


'S'e = {y = y±), 


3,. \ rB_„r,^ . „ . 3^1, 


I J {B,y{A,r - ^Yr) + ^(3A., + B,) + ^ 


]dr - A,y{Y - 4i/el) - ilA^ry 


+ ^ ( e±AA - H{e±,rf + e±,M - Y),r ] - 1 dr 


dr 


AY 


2a|0o 


(2.33) 


Taking the trace of the junction condition and using Eq. (2.13) and the formulas for the gauge transformation, we 
obtain the equation which determines 


(j) 


-A^l= +AAsf, 


(2.34) 


where we introduced 




~ T:^,r ) + —B^yr{3A^r + ^,r) + 


3 J{B,y(^A^r 
+ {B-Y- 2Hil + t±dr)Ail + ilAY^ + J dr 


3?1.AF 


2a^ 


)dr + T ^dr{B + F + nil) 


{^-<Poe±)\ 

a±(^o 


AF 


V - </>o?± 


(2.35) 


Let us consider the junction conditions for the scalar field. Integrating the equations of motion for the scalar field 
across the branes, we obtain the junction conditions for the scalar field as 


±200 — F(£)(0o) • 


±20^^^ = )(0o) 




(0('))V(';y0o) 


(2.36) 


at y = y±. By using Eqs. (2.13), (2.15), and (2.26), the junction conditions for the scalar field in the Newton gauge 
is obtained as 


±200 = F('±)( 0 o), 

a {J) Q 

— - 00 Id = T^^AF^-^) ± e^Asf^ 

A± Ka^0o 


(2.37) 


where we have defined 


A± := 


2 

D±) T 2(0o/ 00 ) 


(2.38) 


and 
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sfun = - ^0 ?V,±, - 5 (»>“' - e)“Ki') ±2 -a,(ivP> - e) - 

+^^Sr + - ^^)] („ = ,±). 

2Ka^(j)Q ^ 2(^q 2«;a^0g2^J 

A more explicit expression is given in Appendix B2. 

Incorporating the boundary conditions (2.32) and (2.37), the master equations become 

[i(TT) ^ J_A](a2A(^)) = 2 k ^ - y^) + , 


(2.39) 


(T=± 


ltd', i]y(d , ^ - »,)A'y> + e<2>Sy . 


(2.40) 

(2.41) 


(T=± 


with 


i^(^)(r) := 2a2 


U) 




3A^Ay( ^) g(./) /A, 

2Ka2(^5 


o 0(7 _ 

2 ^jun 


2k0o ^ 




-I y=y.T 


Solving these equations for and as well as Eq. (2.34) for , and using the gauge transformation Eq. (2.25) 
with the aid of Eq. (2.5), we obtain A^'^\ and C^'^\ which represent the metric perturbations induced on the 

branes. 


III. GRADIENT EXPANSION 


A. Green’s function 


We can write down the formal solution of Eq. (2.40) by means of the Green’s function, 

=‘2kJ 2 f + 2e('^) f d^x T dyGASA , 

( 7 —±. J J V 4 . 


(3.1) 


'v+ 


where the Ga(x, y; x', y') is the Green’s function for the TT part in the static case. The factor 2 in the second term 
reflects the Z 2 symmetry of this brane world model. In the static case, the Green’s function is given by 


GA(x,y;x',y') = - 


(27r)3 


ik(: 


,2K( X —X 


Na{yfa{y' 


E 


Wi{y)wi{y') 

k2 


^Ki 


where Wi (y) is the mode function, and its orthonormal conditions are given by 

ry- rv- 


Widy = 0, 


WiWj 

—^dy = 6^j . 




'y+ 


The normalization factor N is defined by 


N := 


ry- 

2 / a^dy 
Jy+ 


-1 


(3.2) 


(3.3) 


(3.4) 


The first term on the right hand side of Eq. (3.2) is the contribution from the zero mode whose four-dimensional mass 
eigenvalue is zero. The second term corresponds to the propagator due to the Kaluza-Klein (KK) excitations whose 
zth excitation has the discrete mass eigenvalue rriKi- We refer to these modes as KK modes. 

As for the scalar-type perturbations, it has been proved in linear perturbation by considering the source free equation 
that there is no physical mode with a zero eigenvalue of the four-dimensional D’Alembertian [7]. This means that 
the massless scalar-type mode disappears when the stabilization mechanism is taken into account. The explicit mode 
function for the lightest mass mode is found in Ref. [7]. 
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B. Transverse-traceless perturbation 


1. Temporal component 

By the zero mode truncation, in which we substitute only the first term in Eq. (3.2) into Eq. (3.1), we obtain 


A4^Hr,2/) = -2nN ^ atiT^ - + 2Ne^^'> 


<7 — ±. 




^SAdy - 


(7—:t 


(3.5) 


where we did not assume any truncation for the source terms and Sa- We assigned the label 0 to indicate the 
zero mode truncation. 

To evaluate the contribution from the second term in Eq. (3.2) for the TT part of the metric perturbations, we 
follow the strategy that is used in Ref. [7]. Rewrite the part coming from KK modes in the Green’s function as 


d^k 

(^ 


3Zk(3< 


E 


W^Wi 


'‘Ki 


k2 




WiWi f dk^ 


^Ki 


'‘Ki 


{2'kY m 


^ '^‘ki + 


(3.6) 


Under the condition that k'^/"rriYi 1 holds, the first term on the right hand side gives the dominant contribution. 
Notice that the first term is nothing but the Green’s function for Eq. (2.14) with A = 0. Thus to pick up this part of 
the Green’s function is equivalent to solving the equation for by setting A = 0 from the beginning. Substituting 
jYA = where Aq is the zero mode part and is the KK mode part, into Eq. (2.40) and neglecting the 

A term for the KK mode contribution, we obtain 




<7 — ±. 


Nal -S{y- ya) 




where 


ry- 


Sa-=Sa- 2N a^SAdy. 


'y+ 


Applying the integration operator a ^ dya ^, this equation is formally solved as 


4^\r,y)=-2KNj2al^i^^ 


(T=± 


dy'^-Ca +e 


AJ) 


dy' 


dy"a Sa — D(r) 


'y+ 


'v+ 


(3.7) 


(3.8) 


(3.9) 


where (7+ and C- are constants and D{r) is a function of r, which are determined later. The function v±{y) is defined 

by 


v±{y) ^ [ dy'a{y'f . 

a JyT 


(3.10) 


From the orthogonality of Eq. (3.3), must be orthogonal to the zero mode function. This condition fixes the 
constants C± and the function D{r) as 


ry- 


C„=2N 


D{r) = 2N 


dyoA 


, t 'Ja 

dy , 


'y+ 

ry- 


'y+ 


'y-o 


- r dy"a^^A. 


j 2 i 

dya 


(3.11) 




'y+ 


Integrating by parts and taking y = y±^ Eq. (3.9) is reduced to 


(T=± dy+ 

+2NY-^^ r dy'^ r dy"a^SA. 




y+ 


(3.12) 



where we used Eqs. (Bl) and (B3). 

The approximation method that we have used is a kind of derivative expansion method, in which the typical 
wavelength of perturbations is supposed to be long. This expansion is valid only when the smallest mass of the 
KK excitations is sufficiently large. We can obtain the higher order corrections by iteration, in which the A term, 
which we have neglected in the above discussion, is incorporated as the source term. For the TT part, we do not 
consider further iteration than Eq. (3.9). As we will see later, in the scalar-type perturbations, we need to consider 
one iteration of the A term to obtain results accurate to the same order. 


2. Spatial component 

Let us turn to the spatial components of the TT part. Since each spatial component depends on the gauge choice 
of , it is convenient to deal with the gauge invariant combination [38] 


dr{rC^^'>) - 


= —rdr 


+ 2H th + 


{y = y±), 


(3.13) 


where we used Eqs. (2.5) and (2.25), and Sb is given by 






3 • 




Sb = + 2«0o^) -\[B + i{Y + 217^)]- -H{e? + ^ + -5. 


vY 




f dv 


rB, 




V — —pv 


r L 
r\2 


- {C? 


e +(C) - 


(r 


(3.14) 


To fix the gauge degrees of freedom corresponding to the choice of the radial coordinate, we adopt the isotropic gauge 
that is defined by i? = (7 on each brane. With this choice of gauge, the left hand side of Eq. (3.13) becomes rdrB, 
and we immediately obtain the spatial component B (= C) in Gaussian normal coordinates: 

A.B(^)(r,y±) = -iAA(^)T^47’F-^Ar(^)-^e(^) QA5 b + 5^ + 2iLa2±^±j . (3.15) 


( 1 ) 


It is also necessary to specify the explicit form of the radial gauge parameter at first order, , since the second order 


( 1 ) 


perturbations and B^^'> depend on it. Substituting Eq. (2.25) into the isotropic gauge condition, is determined 
as 


( 1 ) 


i± (r) = -^B(^^(r,y±). 


(3.16) 


C. Scalar-type perturbation 

As mentioned earlier, there is no zero mode in the scalar-type perturbation. To evaluate the contribution from 
massive modes, we apply the same technique that we used in the preceding subsection for the KK modes of the TT 
part. First we consider the equation for V'-B by setting A = 0. The homogeneous solutions of Eq. (2.15) with A = 0 
are given by 


u±(y) = 1 - 2B(p)v±(i/), (3.17) 

where v±(y) is defined in Eq. (3.10). The Green’s function Gy, which satisfies 

=-%-y0, (3.18) 

a202 V a2 / 

is constructed as 


Gviyiy') 


3Na\y)a^iy') 

n 


u-{y)u+{y')9{y' 


y) + u-{y')u+{y)e{y - y') 


(3.19) 
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Then using this Green’s function, Eq. (2.41) with A = 0 is solved as 




-N ^ (JUa{y) 


(T=± 


Ki-J) + g(^) 


3u—^Sy 

n^l 


dy' 


(3.20) 


We assign the label 0 because this term gives a contribution to the metric perturbation at the same order as the zero 
mode in the TT part of the perturbation, although it is related to massive scalar-type modes. 

Using Eqs. (B5) and (B17), it becomes 


kN 


AyU(r,9) = -^ i: o.iTy'u,{y) - - 2jvy; 

a=± a 

-2Ne^'^^'^aUcr{y) a^S? -I- [ dy(a'^v-cr/^^cp + - a'^S^) 

a I dy. ^ / 


(3.21) 


where 


4^)(r) := 


HiVa) 


2K<j>Q 


V2^o 


00 

203 


-b 


3(fAY\' 

2Ka20gA 


y=y. 


(3.22) 


As we can see from the first term of Eq. (3.21), this mode of the scalar-type perturbation partly gives long-ranged 
contributions, as we anticipated. We refer to this part of the scalar-type perturbation as the pseudo-long-ranged 
part to distinguish the remaining short-ranged correction. We use the subscript S to represent the short-ranged part, 
although it is also used for the KK mode. 

The source term for the next order correction Yg is given by AIq, which we neglected in the calculation of the 
pseudo-long-ranged part. Since the Green’s function Gy is already known, we easily obtain 


Yy\r,y) 



(3.23) 


Setting y = y± and using Eqs. (B3) and (B6), the expressions for Yq and I 5 are summarized as 

y±) = -^ E T ^ (4 + 2HalSf) - 2Ne^^^ \ ^ aajS'l 


(T—zk 


G=± 


+ r~ dy(a\±AS^ + ^ASy - - 2A , 

Jy+ ^ ‘dn4>l 


wG). ^ d,N fy~ u±AYy\ , 
As (uy±) = - / - ■ 

Ki 


'y+ 


(3.24) 

(3.25) 


D. Large coupling limit 

In the preceding sections, we derived the formal solutions to evaluate the second order perturbations. However, the 
result is very complicated. To simplify the analysis, we assume ^ |0o/0o|) and take the limit 

A± ^ 0. (3.26) 

In the case of the Goldberger-Wise stabilization model [13], this limit corresponds to their large coupling constant. 
In this limit, the junction condition (2.37) for J = 1 becomes 

^W(r,y±)«0oM • (3.27) 

Here we mention that the source term which is given by Eq. (2.37) contains and Hence does 

not vanish even in this limit, and it is reduced to 
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(3.28) 


where we used Eq. (3.27). Therefore the junction condition for the scalar field in this limit is summarized as 

« 0 {y = y±). (3.29) 

Under this condition, the last terms in Eqs. (2.33) and (2.35) vanishes. In particular, Eq. (3.28) gives 

3(pAY \ 

2<^0 y = y,' 

and then we obtain the approximation for (3.22) 

« 0. (3.30) 


IV. RECOVERY OF THE 4D EINSTEIN GRAVITY: FIRST ORDER 


A. Linear perturbation 


We review the results for linear perturbations in terms of the notation of the present paper. From Eq. (3.5), the 
zero mode truncation of the TT part is given by 


a—±L 


B, 


( 1 ) 




(4.1) 


(T=± 


where we have introduced the Newton potential by 

A$±(r) := 47rG p±\r ), 

and G is the induced four-dimensional Newton’s constant defined by 

SttG := kN . 

From Eq. (3.21), the pseudo-long-ranged part in the scalar-type perturbation is obtained as 

(T=± a—±. 


(4.2) 

(4.3) 


(4.4) 


and the gauge transformation is 


e^ = T- 




(4.5) 


In this paper we concentrate on the gravity on one of the Z 2 symmetric branes that carries matter fields on it as 
a source of gravitational field. We assume that the energy-momentum tensor of the matter fields on the other brane 
vanishes. By this simplification, the sum of <I>+ and <I>_ is replaced by 


(4.6) 


in each case. Here we note that, to avoid confusion in showing the formulas for two different situations simultaneously, 
we are using a different convention for the physical length scale from that used in Ref. [7]. 

Substituting the above formulas into Eq. (2.25), we obtain 

4±(’^iy±)=2$±, (4.7) 

where we attached a subscript ± on the perturbation quantities to specify in which case we are working. For instance, 
A+ represents the value of A when only the matter fields on the positive tension brane are taken into account. The 
remaining metric functions turn out to be 

= (7^2 = -2$± (4.8) 

in the isotropic gauge (3.16). These results coincide with the results for the four-dimensional Einstein gravity. 
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B. Correction to the leading term 


To obtain an approximate estimate for the corrections due to the KK mode or the short-ranged part of the scalar- 
type perturbations, it is useful to consider cases in which the back reaction of the bulk scalar field on the background 
geometry is weak; namely, 

1^ ^ 

i/2 3ij2 ’ 

is not as large as unity. In a weak back reaction, the metric is approximately given by the pure anti-de Sitter form 

a{y) « , (4.9) 

and we set 

2/+=0, y-=d. (4.10) 


Here £ is the curvature radius of AdSs. 

If we substitute this warp factor, Eqs. (3.10) and (3.17) are approximated as 

u±{y)-^, v±{y)-U^-l). (4.11) 

2 V / 

Here one remark is in order. The above expression for is not a good approximation near the positive tension brane 
because u+{y+) depends on the difference between N defined in Eq. (3.4) and H y = y+. The value of N in the 
weak back reaction limit is H{a^ — and the difference between H and N is hierarchically suppressed. However, 

unless we consider an extreme case, the deviation of the value of N from this limiting value is not hierarchically small. 
As a result, we have 


u+{y+) = 0(1 ), 


(4.12) 


instead of ©(a^). This also means that the y dependences of and are different. If the single mode with the 
lowest mass eigenvalue dominates the scalar-type perturbation, the y dependence must be the same for both cases. 
Therefore, we find that the modes with higher mass eigenvalues also contribute to the behavior of near the 
positive tension brane. 

Let us consider the KK mode contribution (3.9) in the linear perturbation. A straightforward calculation shows [7] 


^s±(A2/) 


3 



2 



(4.13) 


where we have assumed d/£ ^ 1. Hence, on the brane where the matter fields reside, it becomes 


4+(a?/+)^ A4> 

A^gl{r,y-) « -^A$_ . 


+ ) 


(4.14) 


To compare the KK mode contribution with the zero mode one, we evaluate the ratio between them. Then we find 
that the KK mode contribution is suppressed by the factor 


P± := 


£2 


4 

a^r 


2 

★ 



y+), 


{y = y-), 


(4.15) 


where we introduced a typical length scale r*, and performed a replacement like A « r. 

On the positive tension brane the KK mode contribution is suppressed at r* ^ £ at linear order. Note that, if 
one takes the limit d/£ oo, the KK mode (4.14) seems to diverge. In this limit the lowest KK mode mass goes to 
zero, and the mass spectrum becomes continuous. Then, our expansion scheme which we call a gradient expansion, 
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is no longer a good approximation. Hence, this divergence in the large d limit is just due to the breakdown of our 
expansion scheme. 

On the negative tension brane the KK mode becomes dominant only at the length scale .$0.1 mm when the AdS 
curvature length £ and the hierarchy a+ja- are set to the Planck length £pi and 10^®, respectively. One may think 
that the deviation from four-dimensional Einstein gravity at the sub-millimeter scale provides an observable effect. 
However, in the KK mode contribution to the gravitational potential, <i)± appears only in the form of (dTrG)”^ A$±, 
which is equal to the matter energy density p±. Hence, the KK mode does not contribute to the force outside the 
matter distribution. Therefore, this effect appears to be hard to observe. 

Next, we consider the short-ranged part of the scalar-type perturbation. The short-ranged part (3.23) in the weak 
back reaction limit is evaluated as 


Ys^\r,y) 


dAala+AYg 


( 1 ) 




dy 




-Ayo(r, y), 


(4.16) 


where « const is used (Appendix B3). Here we introduced the lowest mass eigenvalue in the scalar-type 

perturbation ms '■= a-fhs, whose order of magnitude is determined by the last equality [7]. We refer to ms as 
the radion mass. The reason why ms defined above gives the lowest mass eigenvalue can be understood as follows. 
Suppose the mode with the lowest mass squared dominates perturbations in the long-wavelength limit. Then the 
propagator for the scalar perturbation should be proportional to 1/(A — m|). In our approximation of a gradient 
expansion, this massive propagator is expanded as (l/m|) -I- (A/m|) -!-•••. The first term gives and the second 
. Hence, the ratio between them is A/mg. However, again this simple-minded estimate is not correct for Ig^j^ near 
the positive tension brane for the same reason that the approximate expression for u+ [Eq. (4.11)] is not valid near 
the positive tension brane. In fact, the value of Yg^ on the brane can be evaluated by using Eq. (3.25). Substituting 

the estimate given in Eqs. (B24), we obtain Yg^_^{r,yp) « AY^^/rhg. Taking this into account, we guess that the 
formula (4.16) should be modified as 


^i± (a y) 


rfigO^a 


-A$±(r). 


(4.17) 


We give a justification of this formula in Appendix B3. Then, the ratio Kj2(r, ?/±)/Ag2(r, ?/±), where we compare 
the short-ranged part to the Newtonian potential, is 


7± := 


1 


4-2 

a^rrigr 


2 

★ 



(4.18) 


When the radius stabilization mechanism proposed by Goldberger and Wise works most efficiently, the mass ms 
becomes 0{£~^) [7,13]. In this case the short-ranged part of the scalar-type perturbation is suppressed for the same 
reason as the KK mode. We have shown than the zero mode and pseudo-long-ranged part reproduce the correct 
four-dimensional Einstein gravity. The remaining KK mode and the short-ranged part accompany extra suppression 
factors P± and ■j±, respectively. 


V. RECOVERY OF THE 4D EINSTEIN GRAVITY: SECOND ORDER 


We discuss the second order perturbations in the large coupling limit discussed in Sec. HID. As in the case 
of the first order perturbation, we iteratively solve the equations of motion by using gradient expansion. In the 
equations we derive below, we neglect the terms that are relatively suppressed by the factor of 1 /r/ compared to the 
leading contribution. For convenience, we quote the contribution from the zero-type coupling, which is obtained by 
substituting Eqs. (2.12), (2.23), (2.34), (3.5), and (3.24) into Eq. (2.25): 


A4"2(r, y±) = 8^G(p2' + - A ± 2NaiiSf - S^) 


.(2) 


42) 


+2N 


dy 


'y+ 


\Sa± -Sp±)+ a\±A§^± + ^ASy± 

2K0g 


(5.1) 


Here we have used the fact that and Sy: on the brane without matter distribution vanish, which is easily verified 


( 1 ) 

just by noticing that 0 and A±^^(yp) = 0. The contribution from S'-type coupling is 


i(i) 
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(5.2) 


y±) = A y±) - (»’> y±) 


with A^p(r,y±) given by Eq. (3.12), and Yg^\r,y±) by Eq. (3.25). Once we know and Yg^\ the spatial 
component of the metric perturbations is obtained from Eqs. (3.15), (3.24), and (3.25) as 


ABf^(r,?/±) = -^AA''^> +47rG±a±.L± 


1 


i(2) 


( 2 ) 


2 

+3A 


— ASb± 


+ [ dy(a^v±ASip±-\ -^ASy± — 

Jy+ ^ 2 


- 

2 -5^ 


(5.3) 


To identify the order of magnitude of various terms in the second order perturbations, we have to keep track of the 
powers of both r* and a_. Terms with additional inverse powers of r* are basically suppressed for long wavelength 
perturbations. However, as we have seen for the KK mode and the short-ranged part in the analysis of the linear 
perturbations, a complication arises due to the existence of a large nondimensional hierarchy factor 1/a-. Here we 
continue to use the convention a+ = 1. The dependences of the perturbation variables on the warp factor and on r* 
are summarized as 


and 


A, 


( 1 ) 

o± 


w(l) 
-'Oi 5 






a 


a| 


4(1) 


4(1) 

^S- 


Y 


( 1 ) 


S± ’ 


‘Psl 


2 9 9 * 


(5.4) 


(5.5) 


In the following subsection, we first evaluate the terms from the zero-type coupling, classifying them into three parts: 
the part to recover the four-dimensional Einstein gravity, the manifestly suppressed corrections, and the unsuppressed 
corrections. The unsuppressed correction is later shown to be canceled by the contribution from the terms of 5'-type. 
We stress that a weak back reaction is assumed only when we roughly estimate the dependence on r* and a_. 


A. To 


Let us consider Aq given in Eq. (5.1). From Eqs. (2.33) and (2.35) with the large coupling limit (3.29), (5”^ — S^) 
becomes 


aiiSf - Si) = 


3a| J drB^yYr-3 j dr|((.Ar-b a^A,j,(r - -h ^^AF - 2iL(|((.)2 


-ba^ 


-4 


{2(//^dr - Ae){Y + 2He) + BA^y - ^AA + - A),, + r(4^.y + Ae).r 

‘2B yA j- -\- -\- 3-i4) J’ dv . 


(5.6) 


As for S(^, the last two terms in Eq. (2.18) are rewritten as 




3FAr 1 


Ak/)/ 2 




(5.7) 


The expression (5.1) starts with terms of (!l(l/rj); hence we start our discussion with these leading order terms. 
Here, to understand the absence of terms of 0 {r/), we need to notice that 9yA(i) and dyB^A do not have contributions 
from the zero mode, and hence they are 0 {l/r‘/). Let us identify the dependence on the hierarchy a_ of each term, 
concentrating on the case that the matter fields are on the negative tension brane. For this purpose, we can use 
Eq. (5.4). As for dyA^'i and dyB^^\ we use Eq. (5.5) instead because the zero mode contribution exactly vanishes. 
The terms in the second line of Eq. (5.1) possess a y integration. This integration is basically dominated by the 
contribution from the neighborhood of the negative tension brane. One exception is the case in which the integrand 
has the quadratic form of the zero mode contribution of the TT variables (Aq and Bq multiplied by as o^Aq x Aq). 
This integration does not have any inverse powers of a_. The other exception is the case in which the integrand 
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contains the factor of U-/K(j)Q. The formulas for this case are summarized in Eqs. (B24) and we find that only the 
terms with the integrand proportional to give a correction that behaves as 1/a^. The other terms are at 

most Hence, we can pick up the terms with a large power of l/a_ just by looking at the behavior of the 

integrand near the negative tension brane. Then, we find that, among the terms of 0{l/r^), the terms associated with 
a single underline or with double underlines behave as 1/al or l/a?_. Since the usual post-Newtonian correction in 
the four-dimensional Einstein gravity is of 0{a^_/r‘l), we expect that the terms with underlines cancel each other, and 
we show that, in fact, this is the case. The terms with a single underline completely cancel each other. For example, 
the term a?HY^ in Sy is canceled with the last term of Eq. (5.7). Here, it is worth mentioning that the cancellation 
occurs separately within the terms of different types: the terms quadratic in the TT variables, those bilinear in the 
TT variables and the scalar-type variables (Y, (/?, and and those quadratic in the scalar-type variables. The terms 

with double underlines do not vanish completely, but they are, in total, combined to terms of 0{aP_lrl) with the 
aid of Eqs. (B12) and (B13) when we consider the long-ranged part. The contributions from the short-ranged part 
cannot be combined to reduce the power of the warp factor, but they are at most 0{^ j After a straightforward 
calculation, the remaining terms give the usual post-Newtonian term in four-dimensional Einstein gravity (Appendix 
C). This result also applies for the case that the matter field is on the positive tension brane because any term 
of OiXjr^') irrespective of the power of a_ was not discarded in the above computation; namely, we obtain in the 
isotropic gauge (3.16) 


a4±(aj/±) = 87rG(4"^ + wf) - 4$±A<i>± + O (^0 . (5.8) 

Next, we consider the terms of ©(l/rf). Again, we begin by considering the case that the matter fields are on the 
negative tension brane. As we did for the terms of ©(l/r*), we can identify the dependence on hierarchy a_ of these 
terms using the estimates (5.4), (5.5), and (B24). Then, the terms with the highest inverse power of a_ start with 
1/a® , which we refer to as F terms. They are given by 




(5.9) 


,Sy i 


= 2iVA 


J (^^(fidrAifi-SYr^Y'j+Y{AY)^r + 2a^Bs,y(^^'j 


v-' 


- a v±Bs,yYr 


where we have taken into account Eqs. (B12) and (B13). The first term in the curly brackets in comes 

from Sy and the second term from S<^. The remaining terms are at most jat_r'^. The relative amplitude of 
these remaining terms compared to the ordinary post-Newtonian corrections is 0{(3-) or Therefore, only 

the F terms have the possibility of introducing a non-negligible correction. However, we will show in the succeeding 
subsection that the contribution from the E-terms is also completely canceled by that from couplings of the S-type. 

Now we consider the case that the matter fields are on the positive tension brane. In counting the order of each 
term with respect to a_, we will notice that the inverse power of a_ can appear only from a contribution near the 
negative tension brane. Furthermore, from the estimates (5.4) and (5.5), we find that the variables in the first order 
perturbation are at most of ©(l/a/.), and such enhanced variables are associated with the factor l/r*. With this 

notion and the estimate (B24), it will be easy to verify that all the terms quartic in 1/r* in AAq^ are, at most, 
0{£‘^/a^r'^); namely, they are suppressed compared to the ordinary post-Newtonian corrections by the factor of 
Oif3+/a‘i) or 0{j+/a‘i). 

The suppression factors that we encounter at the second order are not as small as those in the linear perturbation, /3_|_ 
and 7 +. This is a natural consequence of our approximation of gradient expansion. Near the negative tension brane, 
the conditions that the scale of the spatial gradient is larger than the typical length scales £ and rrig^, respectively, 
become (/3+/a/_) = {£?<C 1 and ( 7 +/a/_) = (l/m|a/.rj) <C 1. Although the correction seems to become 
large when we consider the case with large l/a_, we think that this is an artifact due to the limitation of the present 
approximation. When we do not have a bulk scalar field, it has been proved that the correction to the four-dimensional 
Einstein gravity in the (l/a_) ^ oo limit stays small [6]. 


B. As and Ys 

In this subsection, we discuss the terms and Yg . The contribution of these terms completely cancels the 
correction due to the F terms. 

From Eq. (3.9), we obtain 
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(5.10) 


where we have performed an integration by parts by using Eqs. (Bl), and also we have used again the fact that Sy: 
on the vacant brane is zero, as well as Eqs. (B2) and (B3). 

First, we concentrate on the case with the matter fields on the negative tension brane. The first term in the square 
brackets in Eq. (5.10) is suppressed by a factor of P' compared with the first term on the right hand side in 
Eq. (5.8), and hence can be neglected. The other terms in are quartic in l/r* or smaller. Hence, we have only 

to study the terms that give a correction of 0{P ja^_r'^). Neglecting the terms higher order in l/r*, the contribution 
from becomes 


-4iV^ 


.(y 


B-^yY-^rdr — ‘SY-A-^y 




dr 


O 


Va4 ’ 


(5.11) 


where we dropped the terms proportional to Y + 2H^y because Eq. (B13) shows that this combination becomes higher 
order in l/a_. As for the terms containing Sa in Eq. (5.10), the contribution of ©(^^/(a® r*)) comes from the terms 
with underlines in Eq. (2.19). For these terms, the integral of oASa- can be performed explicitly as 


a^SA-dy = a^(3a^Y_A_^y + J dra^B_^yY_^r -2 J . 


(5.12) 


Here, note that the contribution from the boundary at y = vanishes. Using Eq. (5.12), we find that the last term 
in Eq. (5.10) cancels the leading order contribution from of Eq. (5.11). Then, the remaining parts in Eq. (5.10) 
give 

AAf_{r,y-) = 2NA, T dj/v_ + [ dra^ B_^yY_^r - 2 [ —+ O Y (5.13) 

Jv+ ^ J J ’’’ ’ Ao 

( 2 ) 

The other correction that we have not considered yet comes from F^j. To evaluate the expression presented in 
Eq. (3.25), first we need to evaluate AFq given in Eq. (3.21). Only the leading terms of 0(l/a/:r*) in AFq are relevant, 
and they are evaluated as 


-AFn 


( 2 ) 


(r,y)« J dr (y_ ,.)2 + ,,)2 - 2a‘^B-^yY-^yr 


+2Na' 


Xau.(y) a(^ 


'Y-ip- KLp_ 


dr 


2a Y- rB_ y - 


24>l 


2r®/^ 


-V-c 




(5.14) 


Note that the terms from S'|’ in Eq. (3.21) cancel the terms obtained by setting y = ya after the y integration in 
Eq. (3.21). Substituting Eq. (5.14) into Eq. (3.25), we obtain 


-AfJ/_^( r, y-) = 3Aa| J dydr^ (^3F_,rAF_ - ‘^ip_drAip_^ + ^(v?-,^)^ - 2a'^B_^y 


- Y_drAY_ 


/ 7 / 

+ / dydr^- 


v+ '■ 


Y-dr^Y- + —(ip-(AY-) r + r)- Y- r(y- r 

' ’ J r ’ \ ’ 


/ ,v 

2Hlp_^. 


+0 


V r/ / 


(5.15) 


From Eq. (B12) the terms inside the second integral turn out to be 0{(A jaAr'^). 
Combining Eqs. (5.9), (5.13), and (5.15), we obtain 
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(5.16) 


i^_+A(4"!-yi!.)) = A|-e^_,, + 4iV J dy^u.J 


+0 


(^)- 

V rf / 


After writing the above expression in terms of U- and v-, we can perform the integration with respect to y by using 
Eq. (Bl). Then, with the aid of Eq. (B3), we find that Eq. (5.16) reduce to terms higher order in 1/r* or those of 
0(^^/alrf). Hence, our conclusion is 


F. + A(a 4 - y^) = O . (5.17) 

The weak back reaction was assumed only for evaluating the order of the residual terms. 

(‘2'\ (“2^ 

In the case with the matter fields on the positive tension brane, the corrections both from Ag_j_ and from are 
suppressed by either {(}+/a?_) or ( 7 +/a^) for the same reason as before. 


C. Spatial components of TT part 


Now the evaluation of in the isotropic gauge is straightforward. Substituting the first order quantities and the 
result of into Eq. (5.3), we basically obtain 


AB(2)(r, y±) = + 4$±A$± - (4>±,r)^ + • ■ 


(5.18) 


which is identical to the result for the four-dimensional Einstein gravity in isotropic coordinates except for the residual 
denoted by (• • •) [38]. 

In the case with the matter field on the negative tension brane, these residual terms are 


Fss- + - Ari?) + , 


where we introduced 


Fsb± = 


-A 


dr 




(5.19) 


(5.20) 


In the same way as for cancellation occurs for the leading order in l/a_ as 

2u- ((/?_,j.)^ 


rv- 


Eq. (5.19) = A i Fsb- + 3N f dydr 

[ Jy+ 


— a V-B-^yY-^. 


O 


atri 


= a 


atrt 


(5.21) 


In the case with the matter fields on the positive tension brane, the residual terms represented by (• • •) in Eq. (5.18) 
are, at most, 0(£^/a?.rf) as before. To conclude, the four-dimensional Einstein gravity is approximately recovered 
under the assumption of the large coupling limit. The corrections to four-dimensional Einstein gravity are suppressed 
by the factor of 0(/3±/at) or 0(j±/a^). 


VI. SUMMARY 

In this paper, we have considered the second order gravitational perturbations in the RS two branes model with 
the radius stabilization mechanism. As a model for the radius stabilization, we have assumed a scalar field that 
has a potential in the bulk and a potential on the brane. From the five-dimensional Einstein equations, the master 
equations for the TT part of the metric perturbations and for the scalar-type perturbation are derived assuming static 
axisymmetric configurations. We have presented formal solutions of these equations by means of the Green’s function. 
We have shown an iterative scheme to obtain approximate solutions by applying the derivative expansion method for 
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the massive modes. For the validity of the derivative expansion, the physical mass of the massive modes must be 
sufficiently large on the respective branes. This sets a limitation on our scheme. Taking infinite separation distance 
between two branes is beyond the framework of the present analysis because the mass of the lowest KK mode becomes 
zero. 

We have shown the recovery of the four-dimensional Einstein gravity in the second order perturbations in the 
following limit: (1) The coupling between the scalar field and the branes is infinitely large [see Eq. (3.26)]. (2) We 
consider the perturbations induced by the matter fields on one brane where we reside, and neglect the effects caused 
by the matter fields on the other brane [see Eq. (4.6)]. 

When we consider the case in which the matter fields are on the negative tension brane, the correction to the 
four-dimensional Einstein gravity appears at the relative order of (!l((a-|-/a_)^(^/r*)^), where I is the AdS curvature 
scale, is the typical length scale of the perturbation, and (a+/a_) is the ratio between the warp factors on the 
positive and the negative tension branes. When this ratio (a+/a_) is 0(10^®), the hierarchy between Planck and 
electro-weak scales can be explained. With this choice of the hierarchy, the correction to the metric in the linear 
perturbation becomes comparable to the usual Newtonian potential when r* < 0.1mm. However, this correction does 
not give a contribution to the force outside the matter distribution. Hence, it seems to be harmless in reproducing the 
predictions of four-dimensional Einstein gravity. We have not confirmed if this feature remains in the second order 
perturbation, but the correction is suppressed by the above factor compared to the usual post-Newtonian correction. 
Hence, the effect due to this correction is almost impossible to detect. 

When we consider the case in which the matter fields are on the positive tension brane, the correction to the 
four-dimensional Einstein gravity in the linear perturbation appears at the relative order of (!l((t'/r*)^), while the 
correction in the second order perturbation is (!l((a_|_/a_)^(f/r*)^) compared to the usual post-Newtonian terms. 
Hence, it seems that the deviation from four-dimensional Einstein gravity appears at a larger scale in the second order 
perturbation. However, this is very likely to be an artifact due to the limitation of our approximation scheme. 

To give a complete proof of the recovery of the four-dimensional Einstein gravity, further extension of the present 
analysis will be necessary. Here we considered the large coupling limit. It will be interesting to evaluate the dependence 
of the correction on the coupling strength. Furthermore, to take into account the contributions from the matter fields 
on the other brane will be interesting. To investigate these issues, a formulation along the line of this paper will be 
promising. Through this second order calculation, we have encountered many miraculous cancellations. This might 
be due to our possibly bad choice of gauge. We would like to defer pursuing a more simplified derivation to a future 
publication, in which we will discuss the unsolved issues mentioned above. 
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APPENDIX A: 4D EINSTEIN GRAVITY 


We have used the isotropic gauge (3.16) to fix the radial gauge coordinates because it is easy to compare with the 
four-dimensional Einstein gravity. However the calculation of the second order perturbation becomes slightly easier 
by taking = 0, although we do not previously know the corresponding four-dimensional Einstein gravity. In this 
appendix we derive the expression for the result of the metric perturbations in four-dimensional Einstein gravity in 
an arbitrary choice of the radial gauge, which corresponds to the various choice of in Eq. (2.25). 

In general, the radial gauge transformation from the isotropic gauge Ajs to an arbitrary gauge A^r is given by 


4(2) _ 4(2) _ Ai ) 


(Al) 


where the generator is related to the quantities at the first order perturbation by the gauge transformation law, 

■2C,;. (A2) 

On the other hand, these bared quantities are also related to the quantities in the Newton gauge as 


1(1) _ o(l) _ 1(1) _ o(l) 
— Ajg Djg 






yl(i)_S(i) 


(A3) 


where is defined by Eq. (3.16). Therefore, we find that CA is simply given by (^ = C ~ Cis- Substituting this 
relation into Eq. (Al), we obtain 
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if = (A4) 

By this equation, the metric perturbation of the four-dimensional Einstein gravity in an arbitrary gauge is deter¬ 
mined. 


APPENDIX B: USEFUL FORMULAS 

In the calculation of the second order perturbations, we often use some relations and results that are easily derived 
from the original definitions and equations, but we have not derived them explicitly. It is convenient to summarize 
such results, and so we devote this appendix to giving the useful relations and formulas. 


1. u± and 


We gives some properties of the functions u± and v± which are defined by Eqs. (3.10) and (3.17). The differentiation 
of these functions with respect to y are 


dy{a^u±) = —2a^Hv± , dy{a^v±) = , dyV± = u± . 


(Bl) 


The last equation is particularly useful to integrate and Sy with respect to y. From the definition, and u_ are 
related as 


1 


u_ — U_|_ = 


2iVa2 


On the branes u± and v± become 


f ^ _ 1 _L Hiy±) 
u±{y±)-l± 

v±{y±) = Tx—TT 7 ) 


u±(2/t) = 1 > 
i’±(2/t) = 0- 


(B2) 


(B3) 


From Eqs. (Bl) and (B3), we obtain 


dy u±vl. = 


1 


'y+ 


2AN^al ■ 


The sums of -I- and — modes are 




H 

7^’ 




^ aua{y)u-aiy) = 0 , ^ aua{y)v-a{y) = 


2Na? 


(B4) 


(B5) 


and also 


XI ^Ay±)fiycr) = ±-^^fiy±) + X ^ 


(T=± 


(T=± 


^ X = Hiy±)f{y±) + Af X ’ 


<7—± 


'^(yv±{ya)f{ya) = - 


f{y±) 


(7 — ±L 


2Na'i ’ 

u±{y^) = Ua{y±). 


(B6) 
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2. Equations 


From Eqs. (2.14) and (2.25) the derivatives of with respect to y are 

‘Ik-- 

= A^y-, (B7) 

_ /\ Ai^,, , 

y) = —{^ - A)- 4:HA^y - —^o{ip - ^oS,y), (B8) 

Qj O 

and using this result we obtain 

A 1 _ _ 2k • •• 

dy{eA,y + = A,y{Y - 4i7e) +e-{Y-A)- + CA,ry - ^(</>oi^ + - 4^) • (B9) 

da o 

We have often evaluated A y and B y on the brane, which are determined by the junction condition (2.32). These 
quantities are purely KK mode contributions. Equation (3.9) gives 


y) = ^ (r, y) = ^J2 > 


3ra^ 


and, taking y = y± , 


4.i(fty±) = , Bg)(r,y±) = , 


(BIO) 


(Bll) 


which are the same as Eq. (2.32). 

As for the scalar mode, we obtain from Eqs. (4.4) and (4.5) 


^0 ^ (a y) + T-V^o (a 2/) = o XI ’ 


<('0 


(>'„<«+2ifa),.,=5E«- 


y=y± 3 


From Eqs. (2.13) and (2.15), 

= 4».(“T.^) + ■ 

3 ' ’ " V / Ina^cpl 


A> = 


3Ay 


2K.a^(po 

Aif 6HAY 


0oV + #^ 


— -1" "I-TT-^ ^(*^0 — H(j)o)Y - 


2 k ; 2 I 4'0 \ 


Integrating by parts, we derive 
ry 


'Vo 


dydyi^o- U—fA [ ^ ^—(7(2/)*^?/’ f dya Syp . 

V 2 k 0 S / Jy^ 


For reference, we give the explicit form of Eq. (2.39) for as 


(B12) 

(B13) 

(B14) 

(B15) 

(B16) 

(B17) 


■Asi, - ^ ^ f ± 7 


200 20g 2Ka^4>l '''' ' 00 l2a^ 

30oAr QHAY Alp {^p - 


A±f0o,,, ,2 , 3ei(Ay)., , A^y + o^HYA , 9(Ay)2 
“ + 2.a20o ^ “ 


2Ka20o 4^20203 


2^0? Ao koP'Ao 

SAT 


+ 


dy 


+{v - ^ 0 ^^) 


2Ka203 


200 


.(A) 

V0o2 

00 

203 , 




{(f - 00?^) 


(B18) 
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3. Goldberger-Wise mechanism 


In the text, we assumed that the y integration containing the factor 1 /(/>q is dominated by the contribution near the 
negative tension brane. To justify this assumption, we discuss the behavior of the bulk scalar field (j)Q. For definiteness, 
we adopt the model proposed by Goldberger and Wise [13], in which the scalar field potentials are 

= (B19) 

where ha ~ (p{ya)- M and ja are the mass and the coupling constant, respectively. The scalar field is solved in the 
weak back reaction approximation as 


My) = (B20) 

where vi = + \/4£“^ + M^, V 2 = and 

B2«</>o(2/+)-e-"i%(y_). (B21) 

According to Ref. [7], the mass squared corresponding to the lowest eigenmode in scalar-type perturbation, m|, is 
estimated by 


N 



dy 


1 


1 

Sa^rrig ' 


(B22) 


Since K(j)Q < for the assumption of the weak back reaction to be consistent, is at most 0{£~‘^a^). The dominant 
contribution to this integral comes from the minimum of {a‘^<Po) ed y = yc = (^log[:/27?2/(j^i7?i)])/4-\/l + 
or y = d when t/c > d. For convenience, we define a quantity of 0(£“^) by m| := absorbing the factor a?.. 

Following the argument given in Ref. [7], we obtain the estimate 


N 



dy 




1 

Salajrhl ’ 


(B23) 


where a and (3 are -|- or —. Here we used the fact that a^u± is a slowly changing function for y > yc in the weak back 
reaction case. From this relation, with the aid of inequalities m_ < aZ^ and < 1, we obtain 



The estimate for Eq. (4.17) can be obtained by approximating Eq. (3.23) as 

YP{r,y)^3NuM) [ ^^^^2 ^ > 2 ^=)- 

Jv+ K(l)^ 


(B24) 


(B25) 


APPENDIX C: EXPLICIT EVALUATION OF THE LEADING ORDER 


To derive the leading order of Eq. (5.8), we give the result of the explicit evaluation of Eq. (5.1). We use Eqs. 
(4.1), (4.4), (4.5) and the junction condition (2.32) that is rewritten in terms of the KK mode as Eq. (Bll). Keeping 
the terms of each term is given as follows: 


-A 


e±Ay+e±Ac + HMM 


= -A 




(Cl) 


21 



4 c± _ 16 f 


±2NalS^ = — dr 




A$j 


2* . ^±$±,r- , $±A$ 4($±,^)2 


r^^±dr — 


Sr 


+ 


3r 


- r A$± 


+ $±..(A$±),, } _ _ 4$^A<i>± - ^r(A$±).. 

±g^|{(4>±.r-)^-264>±A$ + 2 J dr ^($±,r)^ + $±(A$±),, 


(C2) 


±2NaiSf = 


3^^ ’’ mal 


dr 


$±(A$±),^ + -(4>±.r)^ 

r 


_ 4>±.rA$± _ A$± (A$±),, 


6r 


V r® 


4j,3 


64 f , 

+ -/* 


r^$±(ir 


4>±$±,r- 4>±A$± 

j,2 3 j, 


(C3) 




20$±,r 6A$± (A$±), 


) -f( 


164 ^ H 


Na^ 


ff|A*-$±A3.i)+o('i) , (C4) 


2A J dya^S^± ^ \ J 


3$±,rA$± 3(4>±,^)2 $±A$± 17$±4>±,^ 


4 4r 4r 

■3A$± n^±,r 10$± 15 / r2$±(ir 

H 


12r2 

r^$±cir|’ 


54>| 

3r^ 


+(^ ^ ^ (^) 


(C5) 


rv- 

2N / dy'a^v±AS^±= A($^) 

Jy+ 


2 , 

^ H \ 

- 

1± , ) 

9' 

< a^N J 


H 4A /■ 2 3 

^^.^Jdya 


Ol 4 

.n 


(C6) 


, ,u±ASy± 4AA($^) 


3N / dy 

•^y+ '^4>i 


dya^u± + 0(^1 . 


(C7) 


Substituting these results into Eq. (5.1), we obtain the leading term of Eq. (5.8). 
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